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ABSTRACT 

A new model for coevolution of galaxies and supermassive black holes (SMBH) 
is presented that is physically based. The evolutionary track starts with an event 
^ that triggers a significant starburst in the central region of a galaxy. In this 

^ model, the main SMBH growth takes place in post-starburst phase fueled by 

^ recycled gas from inner bulge stars in a self-regulated fashion on a time scale 

^ that is substantially longer than 100 Myrs and at a diminishing Eddington ratio 

^ with time. We argue that the SMBH cannot gorge itself during the starburst 

^ phase, despite the abundant supply of cold gas, because star formation is a 

preferred mode of gas consumption in such an environment than accretion to 
Q the central SMBH. We also show that feedback from star formation is at least 

^ as strong as that from AGN and thus, if star formation is in need of being 

quenched, AGN feedback generally does not play the primary role. The predicted 
^ relation between SMBH mass and bulge mass/velocity dispersion is consistent 

^ with observations. A clear prediction is that early-type galaxy hosts of high 

Eddingtion rate AGNs are expected to be light-blue to green in optical color, 
i gradually evolving to the red sequences with decreasing AGN luminosity. A suite 

^ of falsifiable predictions and implications with respect to relationships between 

>■ various types of galaxies and AGN, and others, are made. For those where 

^ comparisons to extant observations are possible, the model appears to be in 

^ good standing. 

Subject headings: Methods: analytic, ISM: kinematics and dynamics. Galaxies: 
^ interactions. Galaxies: evolution, Supermassive black holes 

X 

;h 1. Introduction 

The tight correlation between galactic center supermassive black hole (SMBH) mass 
(Mbh) and the bulge mass ( Mbg) or velocity dispersion [a) in the nearby universe (e.g.. 



Richstone et al. 1998 Ferrarese & Merritt 2000; Tremaine et al. 2002) strongly suggests 



coevolution of the two classes, at least over the Hubble time. In many semi-analytic calcu- 
lations one of the most adopted assumptions, to put it simply, is that active galactic nuclei 
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(AGN) feedback is able to prevent most of the gas from accreting onto the SMBHs and at the 
same time is able to fix most of the "defects" of galaxy formation models such as the shape of 



the galaxy luminosity function and star formation (SF) history (e.g., Kauffmann & Haehnelt 



2000 Croton et al. 2006 Somerville et al. 2008) with the underlying feedback physics pa- 



rameterized. The substantial success in explaining a variety of observations enjoyed by these 
semi-analytic models is indicative of the relevance of AGN feedback. Calculations of the 
coupled evolution of SMBHs and galaxies using three-dimensional hydrodynamic simula- 
tions deploy thermal energy feedback in regions significantly outside of the Bondi radius 
of the putative SMBH that effectively couples to the surroundings to regulate the SF and 
eventually drive the gas away. These pioneering detailed simulations have provided much 
physical insight and appear to be remarkably successful in accounting for many intricate 
observables, including AGN light curves, Eddington ratio distributions and SMBH-bulge 
relation and its scatter, for certain chosen value of the feedback energy strength (e.g., Di 



Matteo et al. 2005 Hopkins et al. 2006). What is hitherto left open in these calculations is 



the physical origin of the adopted energy feedback. One concern is that the derived SMBH- 
bulge relation depends very sensitively on the adopted energy feedback parameter due to 



the strong radiative cooling (e.g.. Silk & Nusser 2010 Choi & Ostriker 2011). Thus, it is 
prudent to seek underlying physical origins for these successful models and, before that is 
achieved, continue to explore alternative models. 

This paper synthesizes an alternative physical model largely based on known physics. 
Before describing our overall model, we shall first, in §2, examine the plausibility of the 
fundamental claim that AGN feedback is primarily responsible for regulating not only SMBH 
growth but also SF. We argue that scenarios invoking AGN as the primary "blowing machine" 
during the intense starburst phase may logically require significant fine-tuning. We then 
describe the evolutionary path from a starburst to an elliptical galaxy, including the coupled 
evolution of star formation and SMBH growth in the ensuing two sections. 

In §3, we show that growth of SMBH during the starburst phase is limited and consti- 
tutes a small fraction of the overall SMBH consumption. The physical reason is that this 
phase is over-supplied with gas such that only a very small central disc is gravitationally 
stable (Toomre parameter Q > 1) for gas accretion onto the SMBH, while all other regions 
are unstable and more conducive to star formation. Since the SF time scale is much shorter 
than the Salpeter accretion time scale, most of the gas forms into stars. The accreted mass 
during this phase is probably limited to a few percent of the final SMBH mass. In §4, we 
point out that energy or momentum feedback from SF is at least as competitive as that from 
the AGN during the starburst phase. Therefore, SF is largely responsible for blowing most of 
the last patch of gas away to end the starburst phase. In short, during the starburst phase, 
the SMBH does not grow significantly and does not play the leading role in quenching the 
star formation. 
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In §5, we show that most of the growth of the SMBH occurs in the ensuing post- 
starburst period, when the bulge/elhptical galaxy is largely in place and SF enters "passive" 
evolution. The fuel for this primary growth phase is provided by the gas recycled back into 



the interstellar medium (ISM) from aging bulge stars, proposed earlier by Norman & Scoville 



(1988) in the context of a central stellar cluster and stressed recently by Ciotti & Ostriker 
(2007) in the context of elliptical galaxies. It provides a relatively "diffuse" (compared to 
the starburst phase) but steady gas supply that, we show, is ideal for feeding SMBH via 
an accretion disc. Meanwhile, SF is the dominant mode for gas consumption in the outer 
region because the accretion is unstable to fragmentation there, even in this phase. Self- 
regulation is at work for the growth of the SMBH during this period and is provided by 
much more robust (compared to energy feedback) radiation pressure induced momentum. 
The amplitude and slope of the resultant SMBH-bulge relation with this self-regulation is 
consistent with observations. 

In this model, the entire evolution from the onset of starburst, due to a gas-rich merger or 
some significant event that drives a large amount of gas into the central region within a short 
period of time, to becoming a quiescent elliptical galaxy (or a bulge of a future spiral galaxy) 
consists of three distinct periods, as summarized in §5.1 and in Figure [2] (1) "Starburst 
Period": merger of two gas-rich spiral galaxies or some other significant event induces a 
starburst that lasts about 10^ — lO^yrs. The SMBH grows modestly during this period. The 
feedback energy/momentum from the starburst, i.e., supernovae, drives the last patch of 
gas away and helps shut down star formation. (2) "SMBH Prime Period" : several hundred 
million years after the end of the starburst, aging low-to-intermediate mass stars, now in the 
form of red giants and other post-main-sequence states, start to return a substantial fraction 
of their stellar mass to the ISM. The SMBH accretion is mostly supply limited in most of 
this period, except during the first several hundred million years or so, and lasts for order 
of gigayear. Because the rate of gas return from stars diminishes with time, the Eddington 
ratio of the SMBH decreases with time and the SMBH spends most of the time during this 
period at low the Eddington ratio (< 10^'^). The SMBH growth is nearly synchronous with 
star formation from recycled gas during this period. The accompanying star formation rate 
is quite substantial, roughly ~ (5 — 10)(M^,/10"'^^ MQ)(t/lGyr)^^'^ Mq yr~^, where t is time 
in Gyr and M^, is stellar mass of the elliptical galaxy formed during the starburst (at t = 0). 
The duration of this phase depends sensitively on the lower cutoff mass of the initial mass 
function (IMF). (3) "Quiescent Elliptical Galaxy": several gigayears after the end of the 
starburst the elliptical galaxy is now truly red and dead - gas return rate is now negligible so 
both accretion to the central SMBH and residual star formation have ceased. It is possible, 
at least for an elliptical galaxy that is not too massive (i.e., Mtot < 10^^ Mq), that it may 
grow a disk. The feeding of the central SMBH in the bulge of spiral galaxy during this period 
is no longer by aging stars, rather by occasional objects (molecular clouds, stars, etc) that 
happen to be on some plunging orbits due to secular or random events. 
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We present some predictions and implications of this model in §6.2-6.9, followed by 
conclusions in §7. Where comparisons can be made between the predictions of the model 
and observations, they appear to be in good agreement. Some additional predictions could 
provide further tests of the model. 



2. AGN Cannot Regulate Star Formation During Starburst 

While the subsequent sections of quantitative physical analysis are independent of state- 
ments made in this section, we shall argue for the assertion in the title of this section with 
logic, in hopes of being able to provide some conceptual clarity to the role of AGN feed- 
back on star formation during the starburst phase. The starting point of the evolutionary 
sequence is a starburst. It may be triggered by a major merger of two gas-rich galaxies or by 
other significant events that channel a large amount of gas into the central region in a short 
period of time. Consider that an event causes a large amount of gas of mass Mgas to land 
in the central region. Physical processes then operate on the gas to produce a starburst, 
accompanied by some growth of the central SMBH, along with some associated feedback 
from both. Extreme events of this kind may be identified with observed Ultra-Luminous 



InfraRed Galaxies (ULIRGs) (e.g., Sanders et al. 1988) or Sub-Millimeter Galaxies (SMGs) 



e.g., [Chapman et al.|[2005| . Theoretical models (e.g, |Silk fc Rees|[l998| [Hopkins et al.||2006p 



have proposed that feedback from AGN is responsible for the regulation of SF and SMBH 



growth so as to produce the observed Magorrian et al. ( 1998 ) relation where the ratio of the 



final SMBH to bulge stellar mass is Mbh : Mbg ~ 2 : 1000. We shall now re-examine this 
case. 

Consider how the infallen gas may be partitioned. Mass conservation requires Mbh + 
Mbg + Mout = Mgas, where Mout is the amount of gas that is blown away from the bulge. 
Clearly, only a very small fraction of the initial gas Mgas can possibly end up in the central 
SMBH, i.e., /bh = Mbh/ Mgas ^ 1- Let us assume that the reason for a very small /bh 
is because the feedback from the central SMBH prevented its own further growth during 
this phase. Since SMBH masses are observed to span a very wide range, it must be that 
this purported SMBH feedback process that regulates its own growth is galaxy specific, 
i.e., dependent on at least some physical variables characterizing the galaxy. A usual and 
reasonable assumption (which we are not advocating at the moment) for that is that either 
the gravitational potential well of the bulge or of the total halo determines the final SMBH 
mass, in coordination with its feedback. 

Does SMBH feedback dominate that of starburst in terms of regulating both SMBH 
growth and starburst? While we will show later (in §3) that the answer is largely no to 
regulating the starburst at least, we assume that the answer is yes to both for the sake of 
continuing the present thought experiment. The simplified sequence of events then plays out 
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more or less as follows. The central SMBH accretes gas and builds up its feedback strength 
until its mass has reached the observed value, then blows away all the remaining gas and 
both SMBH accretion and SF stop abruptly. What might have happened to SF during all 
this time before the gas is blown away? There are three possible scenarios. Scenario ^1, 
the SMBH accretion is so competitive and quick that most of the gas is blown away by the 
SMBH feedback before much SF has occured. That of course cannot have happened, because 
that would be inconsistent with the observed Mbh — Mbg relation. 

Scenario #2, SF precedes at a pace that is in concert with the SMBH feedback such 
that by the time that Mbh = 0.002(Mgas — Mout), the amount of stars formed is equal to 
Mbg = 0.998(Mgas — Mout); the rest of gas of mass Mout got blown away by the feedback 
from the SMBH. This scenario is designed to match the observed Mbh — Mbg relation. 
What remains undetermined is how large /out = Mout/Mgas is. Is it close to 1 or 0? In the 
case /out ~ 1, because (1 — /out) is a small number, there is no particular preferred value 
for it. The potential well created by the eventual bulge stars would be much shallower than 
the original one already created by the residing gas. In other words, the SMBH only knew 
the potential well of the original gas and it would be rather arbitrary how much stars the 
SMBH decides to allow the bulge to have. If one argues that it is the potential well of the 
total halo mass that matters, the SMBH still did not know how to let SF take place at such 
a rate that we have the very tight observed Mbh — Mbg relation for the bulge region. Thus, 
this case also appears to require much fine tuning. Besides, if (1 — /out) is too small, the 
bulge will be too small compared to what is observed. 

The opposite case with /out <^ 1 is at least substantially more stable, since a large 
fraction of the original gas has formed into stars before the remainder of the gas got blown 
away. In this case the SMBH would "know better" the gravitational potential well eventually 
sustained by bulge stars, because it is not too far from that created by the initial gas. Then, 
how did the SMBH know when to blow away the remaining gas left over from SF and SMBH 
accretion? Should the SMBH blow away the gas when /out = 0.90 (an arbitrarily picked 
number for illustration purpose) or should it wait a bit longer to finally blow away the gas 
when /out = 0.10? It may require more energy or momentum in the former than the latter; 
but that can readily be accommodated by a proportionally increased amount of gas accreted, 
in the vein of feedback from SMBH providing the required feedback energy or momentum. 
Since the amount of gas available before /out = 0.90 is blown away in this hypothetical case is 
capable of growing the SMBH to be 900 more massive than observed and the amount of time 
available (cosmological scale) is much longer than Salpeter time, there is no obvious reason 
why the SMBH cannot grow 10 times (or whatever factor) larger to blow away the gas when 
/out = 0.90 instead of when /out = 0.10. How the SMBH has communicated with the bulge to 
ration the gas consumption would be a mystery. Thus, even in this case with /out ^ 1, taking 
it as a given that the SMBH always stands ready to provide the necessary feedback, having 
SMBH feedback to regulate the overall SF in the bulge such that the ratio of the two matches 
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the observation, again, requires a substantial amount of fine tuning. Nevertheless, since it is 
reasonable to expect that the dependence of the outcome, such as the Mbh — Mbg relation, 
on any proposed feedback processes (including those based on thermal energy deposition 
near the galaxy center) is likely a monotonic function of the adopted feedback strength, it 
should be expected that a solution be found such that the observed Mbh — Mbg relation 
is obtained, for some chosen value of feedback strength, at least for some narrow range in 
Mbg- But, until there is clear physical reason or direct observational evidence to support the 
chosen value of the feedback parameter which the solution sensitively depends on, such an 
approach remains to be refined. We will provide an alternative, significantly less contrived, 
quantitative physical mechanism to circumvent this concern of fine tuning. 



3. Starburst Phase: Modest SMBH Growth and SF Shutdown by Stars 

We have argued in the previous section that AGN feedback cannot logically play the 
leading role in regulating SF, in the sense that while some feedback from the SMBH can 
certainly affect its surrounding gas, there is no particular reason why this could provide a 
quite precise (within a factor of a few) rationing mechanism during the starburst phase so 
as to produce the observed relation between the two. We shall now argue for Scenario ^3: 
during the starburst phase the SF is self-regulated and self-limited, while SMBH growth is 
modest, does not need regulation and does not provide significant feedback to star formation. 

We now give a physical reason for why, even though there is a very large supply of gas 
in the bulge region during the starburst phase, the SMBH growth is modest. We will make 
three simplifying assumptions to present trackable illustration without loss of generality. 
We assume (1) for the regions of interest a geometrically thin Keplerian disc dominated by 
the SMBH gravity (at least at the radii of interest here) is in a steady state, meaning the 



accretion rate (Frank et al. 1992): 



M = 37rz/S ,[1 - (ri„/r)i/2] ' ^ Svrz/E, (1) 



is constant in radius r and time, where is gas mass surface density and u is viscocity; the 
last equality is valid because the radii of interest here are much larger than the radius of the 
inner disc note that it is inevitable to form a disk in the central given the rapid cooling 



and finite angular momentum; (2) we adopt the a-disc viscosity (Shakura & Sunyaev 1973): 

u = acln-\ (2) 
where a is a dimensionless viscosity constant for which magnetorotational instability process 



(Balbus & Hawley 1991) provides a physical and magnitude-wise relevant value; is sound 



speed and f2 is angular velocity (equal to epicyclic frequency for Keplerian disc). The Toomre 
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Q parameter of the gas disc can be obtained from Equations ( Tp ): 



Q 



1 



M 



1/2 



71 



BH 



^3/2^9/4 



(3) 



where G is gravitational constant. The slope of the surface brightness profiles of the inner 
region of the observed powerlaw elliptical galaxies, which are assumed to be the product 
of the starbursts resulting from the gas-rich galaxy mergers, has a value concentrated in 
the range —1.0 to —0.5 (e.g., Faber et al. 1997; Kormendy et al. 2009), reproduced in 



merger simulations (e.g., Hopkins et al. 2009). Presumably the initial gas density profile is 
similar to the final observed stellar density profile in the inner regions. For ease of algebraic 
manipulations, we assume (3) the de Vaucouleur mass surface density profile (with a half- 
mass radius r^) but with the inner region at r < rp = 0.07re modified to be a Mestel disc 
as: ^ 

Sof-l for r<rp, (4) 



S,(r) 



where Eq is the normalizing surface density at some radius ro; we will only be dealing with 
region r < rp; the notional nuclear velocity dispersion of the system without the central 
SMBH at r < rp is related to Sq and rg by 



al = vrGSoro. 



(5) 



Subsequent results do not sensitively depend on the exact slope. The total mass of such a 
hybrid profile is equivalent to a truncated isothermal sphere with a truncation radius of 2re 
and velocity dispersion on galactic scales of ag such that 



1.55a. 
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(6) 



Since the dynamical time, say at Ikpc for a 200km/s bulge being only 5 x lO^yr, is much 
shorter than the Salpeter time, it is appropriate to assume that the gas disc is assembled 
instantaneously with respect to accretion to the SMBH when infalling gas lands on the disc. 
Combining Equations (3 ppp ) we rewrite Q as 



Q 



-1/2-1/2,1/2; 
*0.01 ^0.1 ''E 



r5/4 3 -3/4 
'8 "200'pc 5 



(7) 



where ao.oi = a/O.Ol, e = O.leo.i is the SMBH radiative efficiency, Ie is Eddington ratio. 



Ms = Mbh/10«Mq, (7200 = ^T3/200km/s, 



pc 



r/lpc. The value of a is quite uncertain. 



possibly ranging from 10 ^ to 1 (e.g., Hawley et al. 


1995 


Brandenburg et al. 


1995 


Stone 


et al. 1996 


Armitage 1998; |Gammie 2001| Fleming &; Stone 2003 Fromang & Papaloizou 



2007). Setting Q in Equation ([T]) to unity defines the disc stability radius 



0.22ao.of eol^'/f Mg^/V^ot, PC 
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within which Q > 1 and disc is stable to gravitational fragmentation, and outside which 
Q < 1 and disc is subject to gravitational fragmentation to form stars, supported by both 



simulations (e.g., Gammie 2001; Rice et al. 2003) and circumstantial observational evidence 
of the existence of stellar disc at small Galactic radius (~ O.lpc) (e.g.. Levin & Beloborodov 
2003 Paumard et al.|2006 ). The demarcation value of Q between stability and fragmentation 
does not appear to be qualitatively different even if the disc is under strong illumination (e.g.. 



Johnson & Gammie 2003), as might happen to a nuclear gas disc in the starburst phase. 



The disc mass within rq is 

Mq = 9.8 X lO'a-H'e,'/'lfMl^'a^,l (9) 

This is the accretable mass out of the entire bulge region (note that some of the outer regions 
are more random motion supported). This conclusion reached is in good agreement with 



Goodman (2003), who employs somewhat different assumptions than in this study in that 



he assumes local energy balance, while we impose the observationally inferred inner density 
profile to be self-consistent; the good agreement suggests that this result is quite robust, 
insensitive to assumptions made. Taking cue from our own Galaxy, if we assume that the 
initial SMBH mass of the two merging spiral galaxies of mass ~ 10^^ Mq each is 2.5 x 
10^ M0, for a spiral galaxy of velocity dispersion of 200km/s, we see that the amount of 
mass that could be readily accreted according to Equation ^ using Mg = 0.05 is 6.7 x 
'^J^Mq. Note that the final SMBH mass for such a system is ~ 1.3 x 10^ M0 



10^ao.o(^eo,i 



(Tremaine et al. 2002), if we were to match the observations. 



It is possible that the mass accreted to the SMBH may be larger than that indicated 
by Equation (|9| due to replenishment. Replenishment of low angular momentum gas during 
the starburst phase may be possible in two ways: (1) through orbital decay of outer disc 
gas or (2) direct infall of low-J gas from outer regions. We will show below that (1) does 
not significantly increase the accretable mass. Process (2) is probably unavoidable to some 
extent but unlikely to be frequent enough to be significant for the following reasons. All 
the low angular momentum infalling gas falls into the inner regions initially according to its 
respective specific angular momentum driven by the torque of the trigger event (e.g., merger 
or some other significant torquing event). To replenish low angular momentum gas directly 
to the central region some frequent and significant torquing events during the starburst phase 
are needed. It seems unlikely that such events will be frequent enough to be able to reach 
the observed final SMBH mass: about ~ 100 — 1000 replenishments will be required. One 
might approximately equate the number of replenishment (i.e., significant disturbance) to 
the number of generations of stars formed during the starburst phase (by assuming that each 
generation of star formation manages to redistribute the angular momentum of a significant 
fraction of the gas), which is unlikely to be close to ~ 100 — 1000. In summary, taking into 
account possible additional accretion due to some replenishment and giving the benefit of 
the possibility of ao.oi < 1, it seems improbable that the SMBH is able to acquire a mass 
during the starburst phase that would be much more than 10% of the final value. 
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At r > rg, the disc is unstable to SF. For SF under the conditions relevant here both 
the dynamical and cooling time are short and do not constitute significant bottleneck; if they 
were the only time scale bottleneck, SF would be too efficient. A possible bottleneck for SF 
is the time scale to rid the cloud of the magnetic fiux (assuming the SF clouds are initially 
magnetically sub-critical). The main ionization source in the depth of molecular cloud cores 
is cosmic rays (CR). While the exact ionization rate by CR is unknown for other cosmic 
systems, we have some estimate of that for our own Galaxy, CcR.Gai = (2.6 ± 1.8) x 10"^'^ s~^ 
(e.g., van der Tak & van Dishoeck 2000). If one assumes that the CR ionization rate in 
starburst is 100 times (modeling a typical ULIRG in this case) that of the Galactic value, 
considering that the SF rate in ULIRGs is 100 — 1000 times the Galactic value occuring in a 
more compact region and that the CR in ULIRGs may be advected out via fast galactic winds 
(versus slow diffusion in the Galaxy), one may roughly estimate that the ambipolar diffusion 
time is 7 X lO^yr at a density of n ~ 10^ cm~^ using standard formulas for recombinations 



e.g., McKee & Ostriker 2007). This estimate is, however, uncertain. We will again look 



to direct observations to have a better gauge. Gao & Solomon (2004) show, from HCN 
observations, that ULIRGs and LIRGs convert molecular gas at n > 3 x 10'^ cm~^ at an 
e-folding time scale of tsF ~ 2Myr, consistent with the above rough estimate. It is clear that 
SF time scale is much shorter than the Salpeter time of 4.5 x 10''eo.iyr; in other words, when 
gas is dense and unstable, star formation competes favorably with the SMBH accretion with 
respect to gas consumption. Therefore, most of the gas at r > rg will be depleted by SF. 
When the density profile of the disc at r > rg steepens to be Sc,(r) = TiQ{r /rq)'^^"^, where 
Eg is the gas surface density at ~ rg, the disc at r > rg may become stable again. While 
continued accretion supplied by gas on the outer disc is likely, albeit at a much lower level, 
the mass integral is convergent and most of the mass of this outer disc is at rg given the 
density slope, even if the entire outer disc at this time is accreted. 

Thus, it appears that the amount of gas that is actually accreted by the SMBH during 
the starburst phase is rather limited. This new and perhaps somewhat counter-intuitive 
conclusion is strongly supported by available observations of ULIRGs. This conclusion is 
also opposite to most models that rely on SMBH to provide the necessary feedback to regulate 



star formation (e.g. Silk & Rees 1998 Hopkins et al. 2006). Observational evidence is that 



the SMBHs in ULIRGs and SMGs appear to be significantly smaller (an order of magnitude 
or more) than what the Mbh — Mbg relation would suggest (e.g., [Genzel et alT 1998 Ivison 



et al. 1120001 


Ptak et al. 


2003 


2006 


Alexander et al. 


2008 



ULIRGs should become relatively more important for larger more luminous galaxies (see 
Equation [9]), consistent with observations (e.g., Lutz et al. 1998). Starbursts occuring on 
rotating nuclear disc/rings in ULIRGs are also supported by circumstantial observational 



evidence (e.g., Downes & Solomon 1998). 



The overall conclusion that the SMBH feedback has little effect on the amount of stars 
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formed is in agreement with that of DeBuhr et aL (2010) who investigated the radiation 
pressure-regulated SMBH feedback in the starhurst phase of the merger simulations utilizing 
a sub-grid model for SMBH accretion. One specific common outcome between our calcula- 
tion and their simulation is that most of the gas formed into stars, regardless of the feedback 
strength. A notable difference between our calculation and theirs is that their simulation 
resolution, a gravitational softening length of 47 pc, is significantly larger than tq (Equation 
[s]). As a result, it is possible that their simulations do not resolve small scale that sepa- 
rates stable accretion from unstable, fragmenting disc, which is crucial to our quantitative 
conclusion (note that they use the viscosity parameter a = 0.05 — 0.15 that is larger than 
our fiducial value of 0.01, which would yield a still smaller rg, see Equation [9]) . Thanks 
to that difference, we were able to conclude that, even without considering any feedback 
from the central SMBH, the SMBH during the starburst phase does not grow to anywhere 
close the observed final mass, because star formation can more favorably deplete the gas 
that may otherwise accrete to the SMBH, whereas they find SMBH masses to be too large 
even with substantial feedback (note that they use 10 times L/c radiation pressure force 
assuming multiple scatterings of each converted FIR photon). It seems likely that their 
different conclusion may be due to a much higher accretion rate at their resolution scale, 
which we argue does not reflect the actual accretion onto the SMBH, but rather the disc is 
unstable at that scale and mostly forms stars. As we have noted in the previous paragraph, 
observations indicate that the SMBH masses in the starburst phase appear to be smaller 
than the final values seen in quiescent elliptical galaxies by an order of magnitude, consistent 
with our conclusion. Substantially higher resolution (a factor of ~ 100) simulations may be 
necessary in order to realistically and more accurately simulate the intricate competition 
between accretion and star formation. 



4. A Comparison of Feedback Energetics Between Star Formation and SMBH 

Having shown the unlikelihood of substantially growing SMBH during the starburst 
phase, we now turn to a comparison of the energetics of SMBH and SF to show that, 
feedback from starburst itself should play the leading role in shutting down or quenching 
star formation, i.e., promptly sweeping away the final portion of the gas, where needed. 

To avoid any apparent bias against SMBH or a possibly circular looking argument by the 
assertion that most of the SMBH growth takes place in the post-starburst phase (as we will 
show in §4), we shall for the moment generously assume that the entire SMBH growth occurs 
during the starburst phase, to maximize the energy output from the SMBH, when comparing 
the energetics from the SMBH and the starburst. In Table 1, under the assumptions that 
Mbh '■ Mbg = 2 : 1000, a Salpeter IMF for stars and a radiative efficiency of SMBH 
accretion of 10%, energy output from both SF and SMBH in various forms are listed: (1) 
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total radiation energy, (2) ionizing radiation, (3) X-ray radiation in 2 — lOkeV band, (4) 
mechanical energy, which is supernova explosion energy for SF and broad absorption line 
(BAL) outflow for SMBH, respectively, and (5) radio jets. To obtain energy is ergs per 
Mstar formed, one just needs multiply each coefficient in Table 1 by MgtarC^, where c is speed 



of light. The relevant references are Elvis et al. (1994) and Sazonov et al. (2004) for both 



SBtiiLL) and eBH(2-10keV), RanaUi et al. (2003) for e,(2-10keV), Moe et al. (2009) and Dunn 



et al. (2010) for cbhIBAL) and Allen et al. (2006) for eBH(jet) (if one uses energy seen in 
the most powerful radio jet lobes and assumes that they are produced by the most massive 
SMBH, a comparable value is obtained). The entry for the BAL energy is based on two 
cases and very uncertain, primarily due to lack of strong constraints on the location of the 
BAL and their covering factor. 

It is evident that aside from the energy in the form of radio jets and hard X-rays, SF 
is at least competitive compared to SMBH. Heating due to hard X-rays from SMBH via 
metal line or Compton heating affects only the very central region surrounding the SMBH, 
not over the entire galaxy ( Ciotti fc Ostrike"r||2007 ). Within the physical framework outlined 
here, most of the SMBH growth occurs post-starburst and radio jets occur at a still later 
stage in core elliptical galaxies, energy output (or momentum output derived from it) from 
SF in all relevant forms should dominate over that of SMBH. Our argument that radio jets 
occur at a later stage in galaxy evolution is not at present based on a physical model, but 
on empirical evidence. Observationally, it appears that all significant radio jets are launched 



(Balmaverde & Capetti 


2006 


), with a very 


(e.g., 


Evans et al. 


1999 


Ledlow et al. 


2001 



SO's (e.g., |Veron-Cetty fc Vero^[200l| . But none has been associated with elliptical galaxies 



with an inner powerlaw brightness profile slope. It has been plausibly argued that powerlaw 
elliptical galaxies are produced by gas-rich mergers (we adopt this scenario where a powerlaw 
elhptical galaxy is produced following each major gas-rich merger triggered starburst) (e.g.. 



Faber et al. 1997), whereas core elliptical galaxies are produced later by dry mergers of two 



elliptical galaxies where the fiat core is carved out by the merger of the two SMBHs via 



Table 1. 



# 


Form 


SF 


SMBH 


(1) 


total radiation 


e*(rad) = 7 x 10"^ 


eBH(rad) = 2 x 10"^ 


(2) 


ionizing radiation (> 13.6eV) 


e*(LL) = 1.4 X 10"^ 


eBH(LL) = 3 X 10-5 


(3) 


X-ray (2 - lOkeV) 


e*(2-10keV) = 9 x 10"^ 


eBH(2-10keV) = 5 x 10"^ 


(4) 


mechanical 


e*(SN) = 1 X 10-5 


eBH(BAL) = (0.2 - 2.8) x 10"^ 


(5) 


radio jets 


e=^(jet) = 


eBH(jet) = 4 X 10-5 
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dynamical friction (e.g., Milosavljevic & Merritt 2001). Directly supporting this statement 



is the lack of radio jet in available observations of ULIRGs (e.g., Alexander et al. 2010), in 



agreement with other observations that indicate a significant time-delay between starburst 



and radio activities (e.g., Emonts et al. 2006). An independent, additional argument comes 
from the fact that radio jets are highly collimated and, for the most powerful ones that 
are energetically relevant, they appear to dissipate most of the energy at scales larger than 
that of the bulge region, suggesting that, even if one were to ignore the previous timing 
argument, the efficiency of heating by radio jets for the bulge region is likely low and at 
best non-uniform. Weaker radio feedback, observed almost exclusively in galaxies with an 
atmosphere of hot gas, may be able to steadily provide feedback energy but it is too weak to 
be energetically important. Besides, they appear to only operate in elliptical galaxies with 



hot atmospheres (e.g.. Best et al. 2005). 



The amount of supernova explosion energy that couples to the surrounding medium is 
-E-SN = 1 X 10~^M^,c^, which is exactly equivalent to 5 x 10^'^ Mbhc^ used in the influential 



simulations of Hopkins et al. (2006) with thermal AGN feedback, assuming Mbh : Mbg 



2 : 1000. Because the energy output from supernovae is subject to less cooling than that from 
the AGN, since the former is at larger radii and lower densities than the latter, we expect 
that the amount of energy due to supernovae can at least as effectively as that proposed from 
AGN to drive the gas away. Thus, when most of the gas have formed into stars (i.e., the 
bulge is largely in place after ~ 10'' — lO^yr of starburst), the remaining gas should be blown 
away by collective supernova explosions and the starburst comes to a full stop, reminiscent 



of what is seen in the simulations of Hopkins et al. (2006) with AGN feedback. Detailed 



high-resolution simulations will be necessary, taking into account cooling and other physical 
processes, to ascertain the fraction of gas that is blown away. In short, the bulk of galactic 
winds is likely driven by stellar feedback from the starburst. Galactic winds are observed 
and casual connection between SF rate and wind fluxes has been flrmly established (e.g.. 



Heckman 2001; Weiner et al. 2009), lending strong observational support for the argument. 



5. Post-Starburst: Main Growth of SMBH with Self-Regulation 

The previous section ends when the starburst has swept away the remaining gas and 
ended itself. This section describes what happens next - the post- starburst period, the initial 
period of which is also known as K+A galaxies. 

The newly minted (future) bulge enters its "passive" evolutionary phase, as normally 
referred to. We would like to show that this is when most of the action for SMBH begins, 
fueled by recycled gas from aging low-to-intermediate mass stars. Since two-body relaxation 
time is much longer than the Hubble time, it is safe to assume that the stars formed in the 
inner region during the starburst phase remain roughly in place radially. Angular momentum 
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relaxation may also be ignored for our purpose (e.g., Ranch fc Tremaine||1996 ). However, the 



stellar distribntion in the inner region that initially formed on a disc probably has vertically 
thickened substantially and we will assume that they no longer substantially contribute to 
local gravity on the gas disc (within the thickness of the assumed thin gas disc) subsequently 
formed from returned stellar gas. Because stars in the inner regions are already mostly 
rotationally supported, the shedded gas rains almost "straight down" to land at a location 
that their specific angular momentum allows, to form a disc. Obviously, going out radially, 
the rotational support lessens and star formation may occur in a 3-d fashion. But that does 
not alter our argument about what happened at small radii. The orientation of the disc is 
approximately the same as the previous disc out of which stars in the inner regions were 
formed, since the overall angular momentum distribution of stars has not much changed 
in the absence of any subsequent intrusions. The most important difference of this new 
accretion disc, compared to the disc formed during the starburst phase, is that this new 
disc starts with almost no material and surface density increases with time gradually on the 
timescale of hundreds of megayears to gigayear. 

To have a better gauge how the results obtained depend on the assumed inner density 
slope, instead of assuming a Mestel disc as is done in §3 here we present a more general case 
assuming the inner density profile of the form 

S,(r) = So (^) , (10) 



where n ~ [0.5,1] (e.g., Faber et al. 1997; Kormendy et al. 2009). For this case Equation 



(|8j) is modified, taking into account the gradual change of the gas disc surface density with 
time, to be 



^ , „ , 2/3(3-2n) 




(7r(37r)V2)4/3(3-2n) 

(^J^^^J^^-2/(3-2n)^-2/3(3-2n.)(-^-l/3(3-2n)^5/3(3-2n)^-2/(3-2n)^-2n/(3-2n) ^-^^-^ 

where /rec is the total fractional stellar mass that recycles back to ISM and fg{t) the fraction 
of recycled gas that has returned by time t (out of the fraction free)- The process of SMBH 
accretion in this case goes as follows. The SMBH will accrete all the gas within its Bondi 
radius tb over some period of time, as long as > r^, where rs is defined as 

_GMbh .... 

TB = — ^, (12) 

n 

with cr„ being the velocity dispersion of the inner region of the bulge (r < 20pc or so for 
Mbh = 10® Mq). For the moment we ignore any feedback effect from the SMBH. Since tb 
grows with time and tq decreases with time with increasing fg for r > tb that has been 
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accumulating gas, the condition rq > rs niay be violated at some time t, at which point the 
SMBH is cut off gas supply at its Bondi radius and the SMBH will subsequently grow by 
consuming the final patch of gas on the disc within its Bondi radius. Before the condition 
> tb is reached, the recycled gas that has landed outside (time varying) tb continues to 
accumulate (some of the accumulated gas possibly forms stars). Using Equations ( Ilp2 ) we 
find the turning point rg = is reached when 



n 



fvecfg 



with the disc mass within tq = vb, i.e., SMBH mass, being 



(13) 



3(2 -n) 




~G 



(14) 



From Equation (13) we see that {2 — n)/{2f^cc) > 1 for n = [0.5 — 1]. Thus, we simply correct 



Equation (14) by a factor of 2/i.ec/(2 — n) to finally arrive at 



M 



BH 



3(2 



n) 




frecOi^n 

G 



= 1.9 X 10^(2 - n)2«o.oi/s'eo.i 
with the radius when tq = tb = tbq being: 



a.. 



tbq = 34(2 - n)^ao.oi^i5^eo.i 



200 km/s 



200 km/s 



pc. 



M, 



05 



(15) 



(16) 



Equations ( 15 , 16 ) suggest that the SMBH accreted the recycled gas at r < 20( Mbh/10 M0) pc 



or so for ~ 1; it could be substantially larger for smaller Ie- The reason that the acc- 
retable mass is so much larger during this period than the starburst phase is because the 
accretion disc in this period is replenished continuously at a moderate rate such that it is 
stable within a much large radius than the case of sturburst phase with a much thicker (sur- 



face density-wise) disk. Equation (15) resembles the observed Mbh — c relation (Tremaine 



et al. 2002). We argue the resemblance is deceptive, in a general sense, because it hinges on 
a value of « ~ 0.01 or so and Ie ~ 1. As we mentioned earlier, the currently allowed value 
of a could range from 10~^ to 1 and at the moment we do not know what value nature has 



picked to grow her SMBHs. In light of this situation, using Equation (15) to declare victory 



is premature. However, Equation (15) does suggest that there is enough material and time 



to grow the SMBH to the observed value during the post-starburst phase. This is in stark 
contrast with the starburst phase when there is not enough accretable matter even if one 
pushes the viscosity value to the limit (see Equation [9]). 
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A scenario where the allowed range of viscosity value is limited to one side, i.e., a is 
allowed to have values greater than say 0.01, is much less fine tuned. In this case, some 
self-regulation for the SMBH growth will be necessary. This self-regulation for the SMBH 
growth is indeed achievable during the post-starburst phase, as we will now describe. The 
total amount of radial momentum that radiation pressure of the SMBH may exert on the 
surrounding gas is ec Mbh (this is likely a lower bound, neglecting the possibility of multiple 
scatterings of photons in the optically thick regime). Equating e/3cMBH to /rec Mbg(1— /*)'?^esc 
(that is the momentum of the driven-way gas escaping the galaxy) gives 

1^" e/3c (1 + Aec/.)" I^^'™' ^ ^ 

where A is 

(/e.c/0.15)(l-/,) 

^esc i-i c,\ 



A 



:i + (Aec/0.15)/.)/3r7eo.i 2a' 



where /* is the fraction of recycled gas that subsequently re-formed into stars and fesc is 
escape velocity [for an isothermal sphere truncated at virial radius r„, t>esc(?")/2cr = (1 + 
ln(l + r^/r))^/^ at radius r]; /3 is the fractional solid angle that absorbs the radiation from 
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the SMBH; the term (1 + fescf*) takes into account additional stars added to the bulge 
stellar mass formed from the recycled gas. Of the parameters in Equation [l8l f^sc = 0.15 



is reasonable taking into account that about the half of mass return occuring at early times 
by type II supernovae can escape without additional aid; radiative efficiency of e = O.leo.i 



is consistent with observations (Yu & Tremaine 2002 Marconi et al. 2004); some fraction 



of the recycled gas forming into stars is probably unavoidable, since some gas with column 
density greater than Compton column will slip through radiation pressure (see discussion 
below); also includes the (possibly very large) amount of gas at large radii that would not 
have accreted onto the SMBH in the first place even in the absence of any feedback (e.g., 
molecular clouds on the Galactic disk are not being fed to the Galactic center SMBH in a 
consistent fashion); the factor rj (greater than one) takes into account additional stars that 
are not formed from the starburst event. Overall, considering all these balancing factors, a 
value of A of order unity seems quite plausible. Figure [T] plots the relation between Mbh 
and ctMbg predicted by Equation [17] using A = 1. It is clear that it provides a very good 



fit to the observed data. A similar scaling relation as Equation 17 was derived based on a 



different, radio jet feedback mechanism (Soker & Meiron 2010). 



A similar scenario of linear momentum feedback from AGN radiation pressure has been 



considered by Silk & Nusser (2010) to possibly produce the observed Mbh — Mbg relation 



during the starburst phase but they conclude that the radiation pressure is insufficient by 
an order of magnitude to be able to blow the unwanted gas away. The magnitude of the 
radiation pressure and escape velocity requirement considered here are the same as theirs. 
The difference is that here the amount of gas that need to be regulated in the post-starburst 
phase is nearly a factor of 10 lower and further allowance for star formation from the recycled 
gas make possible that the radiation pressure from the central AGN may be adequate to 
self-regulate the SMBH growth so as not to overgrow it. 



We note that Equation [17] would work without much variation if the gas that is blown 
away is uniformly distributed. The recycled gas is expected to be non-uniform. Even if it were 
uniform initially, thermal instabilities likely make the distribution non-uniform. Given that, 
we elaborate further on Equations ( 17][l8 ) and the physical processes of radiation pressure 
driven winds. Some distinction may be made between about 1/3 of the total solid angle 
where UV and other photons are directly seen from AGN and the other /3 ~ 2/3 of the solid 
angle that has a nearly Compton thick or thicker obscuring screen, most of which probably 
stems from the so-called molecular torus (e.g., [Risaliti et ar[[1999 ). For every AMacc of mass 
accreted, roughly ec/vesc'^Macc = 100eo.i(f esc/300 km/s)^"'^AMacc of mass that rain down by 
aging stars could be driven away by the radiation momentum from the AGN. In the 1/3 
opening solid angle some portion of the radiation pressure driven winds will be accelerated 
to high velocities, perhaps in a fashion similar to what is seen in simulations (e.g 



Kurosawa 



& Proga 2009), observationally manifested as broad emission or absorption lines as well as 
outflows seen in narrow lines (e.g., Crenshaw et al. 2003 Greene et al. 2011). 
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A significant fraction of the material may be accumulated in the remaining /3 ~ 2/3 of 
the solid angle (i.e., Type 2 AGNs), including recycled gas that comes from the other 1/3 
solid angle that is too heavy to be accelerated away "on the fiy" by the radiation pressure. 
In this 2/3 of the solid angle, high velocity winds radially exterior to the molecular torus is 
unlikely given the heaviness (i.e., low opacity) of the molecular torus. We discuss some of the 
physics here. To gain a more quantitative understanding, a look at some observed properties 



of the torus is instructive. Jaffe et al. (2004) measured the radius and height of the molecular 



torus of NGC 1068 to be 1.7pc and 2.1pc, respectively. The mass of the SMBH in NGC 1068 
is (8.3 ± 0.3) X IO^Mq (e.g., [Marconi fc Hunt||2003[). If we extrapolate to a 10* Mq SMBH 



assuming that the location and height of the molecular torus is proportional to the SMBH 
mass, we have a surface area of the torus equal to 3200 pc^ at a SMBH-centric radius of 
20pc. If we assume that the column density of the molecular torus is lO^^cm"^ (e.g., Risaliti 



et al. 1999), its total mass is then 2 x 10'' M0. The dynamical time at 20 pc is 10^ yrs. 
A SMBH of mass 10* Mq accreting at Eddington rate would grow a mass of ~ 10^ Mq in 
10^ yrs, while the overall rate of gas return would be ~ 2 x 10^ M© over the entire bulge 
during that period. Thus, the abundant gas supply rate suggests that the necessary (not 
sufficient) condition for a near "steady" state is met such that the molecular torus may be 
kept roughly invariant with time, with the rate of driven-away gas by radiation pressure plus 
that of gas forming into stars equal to the rate of gas return from aging stars. 

Given the short star-formation timescale of the very dense gas in the molecular torus, 
it would be unavoidable that star formation should occur there (as well as some regions 
exterior to it). This "lightens up" the torus to the extent that it may be pushed away by 
the radiation pressure, when the condition that the deposited radiation momentum divided 
by the accumulated mass exceeds the escape velocity (assuming, in the absence of radiation 
pressure, the torus would just be in a bound circular orbit). In this sense the radiation 
momentum from the SMBH serves to retard gas supply to accretion from the torus to let 
SF take over to have it mostly depleted. In combination with the analysis in the preceding 
paragraph, it seems physically plausible that radiation pressure and depletion of gas by star 
formation is able to jointly reduce and regulate the amount of gas that feeds the central 
SMBH. Given that the overall margin, in an "on average" sense, is quite thin (i.e., A ~ 1 in 



Equation 18), it is likely that there are significant variations in A, perhaps up to a factor of 
a few. 



In the 1-d simulations of Ciotti & Ostriker (2007) for an elliptical galaxy, the SMBH 



growth appear to be intermittent. The intermittency in their simulations was caused by a 
hot X-ray heated bubble that prevents continued gas accretion, until it bursts, which is then 
followed by another accretion episode, and so on. We suggest that Ray leigh- Taylor instability 
on the shell enclosing the X-ray bubble may prevent the X-ray bubble from infiating, as 
hinted by recent 2-d simulations of Novak et al. (2010). It it reasonable to assume that 



shell fragmentation in three-dimension is still more pronounced to allow continued defiation 
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of a notional X-ray bubble. Observationally, the lack of significant X-ray emission from 
circumnuclear region in powerlaw elliptical galaxies host AGNs, which we argue are the 
post-starburst galaxies we consider here, supports the picture that the hot bubble is not 
robust (e.g., Pellegrini 2005). In the absence of a hot X-ray bubble guarding the SMBH, we 
suggest that the recycled gas from aging stars is able to reach the disc and the accretion, 
with self-regulation argued above, is quasi-steady without major flares of magnitude seen in 
1-d simulations. As we will show later, a steady declining accretion rate proportional to the 
gas return rate provides a much better match to at least two observations: (1) the observed 
early-type host galaxies of AGNs are mostly in the green valley of the galaxy color-luminosity 
diagram with a small fraction in the red sequence (§5.2) (e.g., Salim et al. 2007 Silverman 



eraL][2008| [Hickox et al^[2009j [Schawinski eFaLjlMol , but very few in the blue cloud, which 



would have been the case if AGN flares are accompanied by starbursts (Ciotti & Ostriker 



2007); (2) the observed AGN accretion rate for early-type galaxies in the local universe 



displays a powerlaw distribution with the amplitude and decay rate (Kauffmann & Heckman 



2009) that is expected from the non-fiare scenario that is proposed here. This indicates that 



bursty AGN accretion, while quite possible and sometimes perhaps unavoidable, is probably 
not the dominant mode. It is currently a challenge but will be of great value to carry out 
3-d high-resolution simulations to more accurately quantify this outcome. 



6. Model Predictions and Discussion 

We have presented a physically motivated picture for the coevolution of galaxies and 
SMBH starting with a triggered starburst. Let us now summarize the entire evolution in §5.1 
and then give an incomplete list of implications and predictions in §6.2-6.9 to be qualitatively 
compared/verified with observations. 



6.1. Three Distinct Periods of Coevolution of Galaxies and SMBH 



From the onset of a significant central starburst to becoming a quiescent bulge there are 
three distinct periods, as summarized in Figure[2]for an example merger of two gas-rich spirals 
each of mass ~ 10^^ M0 that eventually becomes a powerlaw elliptical galaxy of velocity 
dispersion of 200km/s. We stress that the trigger event is not limited to major mergers. 
This three-stage scenario is not new and its successes with respect many observations have 
been discussed previously (e.g.. 



Granato et al. 


2004 


2006 


Cirasuolo et al. 


2005 


Lapi et al. 



2006 Lamastra et al. 2010). The new theoretical element here is the primary growth of 
SMBH in the post starburst phase, which is reflected in the color and other properies of 
AGN hosts and we will show is in remarkable accord with latest observations, in contrast to 
the conventional scenario where SMBH growth primarily occurs during the starburst phase. 
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The time boundaries between difference consecutive phases (three ovals) are approximate 
(uncertain to a factor of at least a few). Given the complexity and variety of starburst trigger 
events, one should expect significant variations from case to case. The expected consequences 
or predictions of this model are in many ways different from and often opposite to those of 



models that invoke AGN feedback to shut down both starburst and AGN activities (e.g, Silk 



& Rees 1998 Hopkins et al. 2006). A new and in some way perhaps the most fundamental 



finding of this work is that the SMBH does not grow during the starburst phase as much as 
previously thought, required in AGN-feedback based models, despite the obvious condition 
that there is a lot of gas being "jammed" into the central region; this is different from almost 
all previous work (e.g. Silk fc Rees||1998 Hopkins et al.||2006 DeBuhr et al.||2010 ) that either 
need to advocate very strong SMBH feedback or appear to overgrow the SMBH. 

The idea of feeding the SMBH with recycled stellar material in the post-starburst phase 
is not new (e.g., Norman & Scoville 1988; Ciotti & Ostriker 2007) and we inherit most of 



the already known elements from prior work, including gas return rate and the likelihood of 
continued star formation. Our analysis shows the likelihood that the SMBH may be fed too 
much in the post-starburst period in the absence of feedback from the SMBH, in dramatic 
contrast with the starburst phase when SMBH feedback is insufficient. While energy feedback 
from the SMBH certainly plays a role, we show that the more robust momentum feedback 
from SMBH radiation pressure can play a critical role in regulating SMBH growth, not 
necessarily only by blowing powerful winds, but rather, in combination, by also pushing 
away thus retarding accretion of unwanted (by SMBH) gas to be instead consumed by star 
formation. While our analysis may have captured some of the essential physics in terms of 
accretion and star formation demarcation, to more realistically model the complex accretion 
and star formation dynamics, much higher resolution 3-d radiation hydrodynamic simulations 
will be required and will be of tremendous value. 

The "size" of the starburst depends on the "size" of the triggering event, with at least 
some fraction of ULIRGs and SMGs due to major mergers of massive gas-rich gas. However, 
irrespective of the size of the starburst event, the time scales involved, being largely due to 
physics of stellar interior and accretion time scale, remain the same. 

(1) "Starburst Period": this phase is triggered by some event. The SMBH grows modestly 
during this period to possibly attain a mass that is up to order ten percent of its final mass. 
This phase lasts about 10'' — lO^yrs for typical starbursts and the host galaxies during this 
phase are in the blue cloud in the luminosity-color diagram. The feedback energy/momentum 
from the starburst, i.e., supernovae, drives the last patch of gas away and shuts down star 
formation, if needed. In other words, the starburst is self-regulated, not by the central AGN 
during this period. 

(2) "SMBH Prime Period" : several hundred million years after the end of the starburst, aging 
low-to-intermediate mass stars, now in their post-main-sequence phases, start to return a 
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substantial fraction of their stellar mass to the ISM. The SMBH accretion is fueled by this 
recycled gas lasting for order of gigayear. The growth of SMBH is self-regulated, readily 
provided by the radiation pressure from the AGN. The host galaxies during this period start 
out light-blue or in the "green valley" and migrate to the "red sequence" . Because the rate 
of gas return from stars diminishes with time and SMBH mass grows, the Eddington ratio 
of the SMBH decreases with time. The SMBH growth is synchronous with star formation 
from recycled gas during this period. The accompanying star formation rate may also be 
substantial but typically does not constitute a starburst during this period. The entire 
duration of this phase depends sensitively on the lower cutoff mass of the initial mass function 
(IMF) - a sensitive and powerful prediction of this model. 

(3) "Quiescent Bulge" : several gigayears after the end of the starburst the bulge is now truly 
red and dead - gas return rate is now negligible so both accretion to the central SMBH and 
residual star formation have ceased. It is possible that a disk is grown later around the 
bulge. The feeding of the central SMBH in the bulge of spiral galaxy during this period is no 
longer by overhead material from aging stars, rather by occasional objects that happen to 
be on some plunging orbits to be disrupted by the SMBH and form a short-lived accretion 
disc. Candidate objects may include molecular clouds, some tidally disruptable stars or gas 
streams. Significant disturbances or torques, such as minor mergers and galactic bars, could 
provide the necessary drivers for some more consistent accretion events. How is a red and 
dead bulge with a hot atmosphere able to remain star-formation-free? This is a major topic 
on its own right and beyond the scope of the current paper, but will be addressed in a future 
paper. 



6.2. Some "Obvious" Implications of the Model 



There are some unambiguous discriminating signatures of this model that already can 
be directly "read off' Figure |2} We highlight several here. 

(1) Starburst and AGN growth are not coeval in this model. AGN does not regulate the 



starburst, consistent with observations (e.g., Schawinski et al. 2009; Kaviraj 2009). AGN 



activities is expected to outlive the starburst, in agreement with observations (e.g., Geor 



gakakis et al. 2008). These predictions are opposite to those of models that invoke AGN 



feedback as the primary regulating agent. 

(2) The apparent requirement of a rapid migration of early-type galaxies from the blue 
cloud to the red sequence, in order to produce a bimodal distribution in color (e.g., Blanton 



et al. 2003), is primarily due to the prompt shutdown of SF by stars (i.e., supernovae) 



at the end of the starburst phase; there is no need to invoke other ingredients, consistent 



with observations (e.g., Kaviraj et al. 2010). Observationally, there is no evidence that the 



- 21 - 



presence of an AGN is related to quenching of star formation or the color transformation of 



galaxies (e.g., Aird et al. 2012). This prediction is different from that of models that invoke 



AGN feedback to quench star formation. 

(3) AGN activities in ongoing starburst galaxies, i.e., buried AGN activities, are not 



expected to be dominant in this model, in agreement with observations (e.g., Genzel et al. 



1998 Ivison et al. 


200C 


); 


Ptak et al. 


2003 


1 


vison et al. 


2004 




Alexander et al. 


2005a|b 


Schweitzer et al. 


2006 


] 


^awakatu et al. 


2006 


Alexander et al. 2008; Veilleux et al. 


2009) 



Note that the above statement is not inconsistent with AGN/QSOs being associated with 
galaxies in the process of merging, which may enhance accretion activities in the involved 



(yet to merge) galaxies (e.g, Bahcall et al. 1997 Hennawi et al. 2010; Smith et al. 2010). 



(4) The most luminous quasars that accrete with high Eddington ratios occur order 
of lOOMyr after the end of the starburst. They may contain substantially more merger 



signatures, which appears to be indicated by observations (e.g., Bennert et al. 2008). If one 
were to identify a population in-between ULIRGs and more regular QSO hosts in terms of 
spectral properties, they should show some more signs of tidal interactions that are yet to 
fully settle since the starburst, also consistent with observations (e.g., Canalizo & Stockton 



2001). 



(5) Low Eddington ratio AGNs that are expected to last order of Gyr are not expected 
to show a close linkage to major disturbances that trigger the starburst (e.g., mergers), since 
possible signatures of the trigger merger event have largely been erased over time, consistent 
with observations (e.g., Grogin et al. 2005 Cisternas et al. 2011). Thus, one does not 



expect to see merger signatures to be associated with moderate-luminosity AGNs, which is 
in contrast with AGN feedback based models where most of the moderate luminosity AGNs 
are expected to coincide with starburst. 

(6) While the green- valley morphologically early-type galaxies that host AGN is the evo- 
lutionary link between starburst galaxies (in the blue cloud) and the red elliptical galaxies 
(on the red sequence), it is useful to distinguish between them and the other class of green 
galaxies that simply continuously form a modest amount of stars (such as our own Galaxy). 
The former are chronologically immediate successors to starburst galaxies and should be in 



early-type galaxies, strongly supported by observations (e.g., Salim et al. 2007 Silverman 



et al. 2008; Hickox et al. 2009 Schawinski et al. 2010), whereas the latter are not a chrono 



logically intermediate class between the blue cloud and the red sequence. The total green 
galaxy population will be the sum of these two different morphological types, with some 
obvious implications, such as green galaxies having mixed morphological types with limited 
merger signatures, consistent with observations (e.g., Mendez et al. 2011). This prediction 



is in contrast with AGN feedback based models where most AGN hosts are expected to co- 
incide with starburst and a small fraction, mostly the most luminous AGNs (occuring near 
the end of the starburst phase), is expected to have matured early- type morphologies. 
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(7) While the early-type AGN host galaxies may have similar morphologies as and will 
eventually evolve to inactive elliptical galaxies, the former should have much bluer colors than 



the latter, consistent with observations (e.g., Sanchez et al. 2004). The basic morphological 
properties of the host galaxies of the most luminous quasars, corresponding to the most 
massive SMBHs in the prime growth phase should resemble those of giant elliptical galaxies. 



consistent with observations (e.g., Dunlop et al. 2003). 



(8) Because of the expected rate of gas return (oc t^^'^ on gigayear scales) to which 
both SMBH accretion and star formation are proportional and because more powerful AGN 
accretion occurs closer in time to the preceding starburst, it is expected that more powerful 
AGNs are hosted by early-type galaxies with younger mean stellar ages, consistent with 



observations (e.g., Kauffmann et al. 2003 Jahnke et al. 2004). 



(9) The accompanying star formation rate of elliptical galaxies may be quite substantial, 
on the order of ~ (5 - 10)(M,/10" M0)(t/lGyr)-i-3 Mq yr'^. Thus, while most AGN host 
galaxies have left the blue cloud, a significant fraction of them, especially those hosting lumi- 



nous AGNs, should still have substantial SFR, consistent with observations (e.g., Silverman 



et al. 2009 Shi et al. 2009). It is expected that the incidence of star formation signatures 



(e.g., dust) in the nuclear region should correlate positively with AGN activities for elliptical 
galaxies, because the strengths of both are proportional to the gas return rate, consistent 
with observations (e.g., Simoes Lopes et al. 2007). These predictions are opposite to AGN 



feedback based models where star formation is expected to be completely quenched after 
AGN feedback clears the gas out. 



6.3. Origin of Two AGN Accretion Regimes 



Kauffmann & Heckman (2009) presented an insightful observational result of two dis- 
tinct regimes of black hole growth in nearby galaxies along with its apparent implications. 
They find that star-forming galaxies display a lognormal distribution of Eddington ratios; 
their interpretation is that in this regime accretion on to the SMBH is not limited by the 
supply of gas but by feedback processes that are intrinsic to the SMBH itself. Our model 
provides the following alternative interpretation for this phenomenon: this lognormal dis- 
tribution merely reflects two random processes at work: (1) the amount of gas that landed 
on the stable accretion disc to provide accretion to the SMBH during the starburst phase 
depends on many "random" variables of the triggering event (in the case of a merger, such 
as merging orbit inclination, velocity, spin alignment, etc), and (2) observations catch a 
random moment during the accretion of this gas. Central theorem should then give rise to 
a lognormal distribution. Another class of possible triggering events for SMBH accretion in 
star-forming galaxies (e.g., dormant SMBH in the bulge of disk galaxies) is stochastic feeding 
due to some random events, which should also follow a lognormal distribution. 
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Separately, they find that galaxies with old stellar populations is characterized by a 
power-law distribution function of Eddington ratios and the AGN accretion rate is about 
0.3 — 1% of the gas return gas from recycling. In our model the expect accretion rate 
is expected to be MBH/(/rec Mbg) = 1-3 x 10~^Ae(^}(j2oo- This expected relation between 
SMBH accretion rate and gas return rate is remarkably close to their observed value. As 



Kauffmann &; Heckman (2009) already pointed out, the powerlaw distribution is consistent 
with the recycling gas return rate oc t~^'^ (Mathews 1989). This is a strong support for the 
proposed model here. 



6.4. Initial Mass Function and AGN Accretion History 



Because the least massive stars live the longest, the cutoff mass of the initial stellar 
mass function (IMF) plays an important role in shaping the evolution on longer time scales 
of > IGyr. For example, a 0.92 Mq star (solar metallicity) has a lifetime of lOGyr, whereas 
a 1.4 M© only lives ~ 2Gyr. Thus, the duration of the "SMBH Prime Period" depends 
sensitively on the lower mass cutoff of the IMF. Figure [3] shows several cases of the evolution 
of the SMBH growth tracks. It shows that the evolution and duration of SMBH growth in 
the post-starburst phase depend sensitively on the low mass cutoff of the IMF. We see that 
for a cutoff mass of 0.92 Mq the SMBH spends about lOOMyr accreting at Eddington limit 
when its mass is up to about 10% of its final mass and a significant period (> IGyr) at less 
than 1% of the Eddington rate, and most of the time at about 0.1% of the Eddington rate 
when its mass approaches its final mass. On the other hand, with a mass cutoff of 1.4 Mq 
the entire SMBH accretion shortens to 2Gyr and does not extend below 10^^ Eddington 
rate. Since not all elliptical galaxies at present time are observed to accrete at 0.1% of the 
Eddington rate, this already suggests that a higher than 0.92 Mq cutoff mass in the IMF may 
be required. Presently there is circumstantial evidence for massive star formation in galactic 



centers, including our own Galaxy (e.g., Lu et al. 2009) and M31 (e.g.. Bender et al. 2005). 



Given the very sensitive dependence of stellar lifetime on stellar mass, careful considerations 
along this line may prove to be very powerful in placing constraints on the low-mass cutoff in 
the IMF as well as testing this model. Detailed comparisons between theoretical prediction 



with observational data in terms of the AGN luminosity-mass plane (e.g., Steinhardt & Elvis 



2010), the Eddington ratio range (e.g.. Woo & Urry 2002), AGN ages at different redshifts 
(e.g., Martini|2004 ) or at different luminosities (e.g., Adelberger et al.|2005 ) should also prove 
very powerful in constraining the IMF. We note that our assumption used to derive the light 
curves in Figure [3] is extremely simplistic and therefore we do not expect that they provide 
satisfactory matches to observations. It is called for that additional ingredients be included 
to account for, e.g., variations in stellar distribution, possible variations of IMF as a function 
of local star formation conditions, dependence of initial seed SMBH mass on galaxy model, 
etc, in order to have a more encompassing analysis. We shall carry out a more detailed 



analysis with additional parameters in a future study, especially when measurements of both 
SMBH masses and accretion rates become significantly more precise for a large sample of 
active galaxies. 



6.5. Super-Solar Metallicity of Accreting Gas 



One clear implication is that the accretion gas, being shedded from aging stars, should be 
very metal rich with supersolar metallicity, in agreement with observations (e.g. Hamann & 



Ferland 1993), especially to explain super-solar N/He ratio (e.g., Hamann & Ferland 1999). 



This is because nitrogen is believed to be secondary nature, where its abundanace scales 
quadratically with metallicity. The recycled gas that is feeding the SMBH in our model 
fits the bill most naturally. In addition, the metallicity of accretion gas is not expected to 
depend on redshift, being intrinsic to stellar evolution, consistent with all accreting gas being 



very metal rich at all redshifts, including the highest redshift SDSS quasars (e.g.. Fan et al. 
2006|. 



6.6. Relative Cosmic Evolution Between Starburst Galaxies and AGN 



Given the modest amount of time delay (several lOOMyrs) between the starburst phase 
and the SMBH prime growth phase, it is unsurprising that one should expect to see nearly 
synchronous evolution between the starburst and SMBH growth on longer, cosmic time 



scales, consistent with observations (e.g., Boyle et al. 1988 Nandra et al. 2005). 



In the context of the observed cosmic downsizing phenomenon, the downsizing of galax- 
ies (e.g., Cowie et al. 1996 Treu et al. 2005) should thus be closely followed by downsizing 



of AGNs (e.g., Barger et al. 2005 Hasinger et al. 2005). There is, however, a very impor- 
tant difference between the two classes in post peak activities, predicted in this model. For 
starburst the shutdown time scale is expected to be about ~ lOOMyrs, whereas for moderate- 
luminosity AGNs (i.e., Eddington ratio ~ 10~^) the decay time scale is of order of ~ IGyrs. 
With a deep AGN survey that is capable of subdividing early-type galaxies in terms of their 
masses, one should be able to differentiate between the downturn time of starburst galaxies 
and that of AGNs hosted by elliptical galaxies at a fixed mass. This prediction would be a 
strong differentiator between this model and AGN-based feedback models. 
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6.7. AGN Broad Emission and Absorption Lines 

Some of the overhead material raining down onto SMBH accretion disc from recycled gas 
from aging low-to-intermediate mass stars provides the material observed as broad emission 
lines (BEL) and broad absorption lines (BAL). When some of this gas, probably in the form 
of some discrete clouds, reaches the inner region of the the SMBH (at r < lOV^, where 
is Schwarzschild radius), the clouds will be accelerated by radiation pressure, likely through 
some absorption lines, to velocities up to 0.1c. These clouds will be the observed BEL and 
BAL. The fact that only 15-20% of type I AGN to have BAL may be indicative of the 
discrete nature of the clouds, not unexpected from discrete stellar remnants or from cooling 
instabilities. 

An advantage of this overhead material is that it naturally provides gas clouds that are 
presumably to be some > 50° off the equatorial plane, in order not to be obscured by the 
molecular torus (there are of course BEL and BAL gas clouds at smaller angles but they 
are not seen directly). In this model we do not need any additional pressure force to lift 
the gas off the accretion disc - some of the raining down gas clouds from aging stars will be 



launched outwards before they reach the disc, physics of which is well known (e.g., Murray 
et al.||1995l). 



6.8. Evolution of SMBH Mass Relative to Bulge Mass 



Massive elliptical galaxies appear to have increased their masses by 30 — 100% in the 
last 7Gyr (e.g.. Brown et al. 2008). The growth of the elliptical mass is not expected to be 
always accompanied by corresponding growth in the mass of the central SMBH. For example, 
merger of a spiral galaxy without a significant SMBH and an elliptical galaxy would make 
the final SMBH appear less massive. Given the dependence of Mbh/ Mbg oc a oc (1 + z)^/^ 
predicted in this model, we predict that the Mbh/ Mbg relation should evolve with redshift 
stronger than (1 -|- z)^/^ for quiescent elliptical galaxies. 



6.9. On Relation between SMBHs and Pseudo-bulges 



It is useful to add a note on the difference between classic bulges and pseudo-bulges 



(Kormendy & Kennicutt 2004) with respect to the central SMBHs in this model. The 



relation derived. Equations (17, 18), that matches the observed Mbh — Mbg relation is 



dependent on the abundant supply of recycled gas in the inner region. Given the sufficient 
gas supply from recycled gas, the feedback from the SMBH then can regulate its own growth. 
This essential ingredient of sufficient gas supply is consistent with the observed inner slope 
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of classic bulges (e.g., Faber et al. 1997; Kormendy et al. 2009), as we have shown. 



The situation would be very different, if star formation is not as centrally concentrated 



as in classic bulges, for example, in rings (Kormendy & Kennicutt 2004, and references 



therein) of high angular momentum with a hollow core. In this case, the amount of recycled 
gas raining down from the innermost region may depend on other unknown factors. For 
instance, if secular processes act promptly, compared to the time scales of stellar gas recycle 
(~ 0.1 — 1 Gyr), to be able to substantially fill the central region with stars initially formed 
in outer regions, the SMBH may follow the track we described. If, on ther other hand, 
secular processes evolve on longer time scales, the recycled stellar gas would predominantly 
land in outer regions that do not efficiently accrete to the SMBH, which would in turn not 
grow substantially. It would seem likely that there may be two trends for pseudo-bulges: (1) 
there will be large variations in Mbh — Mbg relation and (2) SMBH masses may lie below 
that of the Mbh — Mbg relation derived from inactive classic elliptical galaxies/bulges, both 
consistent with independent considerations in the context of hierarchical structure formation 



model (e.g., Shankar et al. 2012). Observations, while very challenging, may have already 



provided some hints of both (Greene et al. 2008). 



Moreover, we do not expect any discernible correlation between the SMBH and galaxy 
disk or dark matter halo, simply because the stars in disks do not affect SMBH growth and 
the overall dark matter halo, while indirectly affect the escape velocity that enters Equation 



(18), does not control the amount of gas that feeds the SMBH. This prediction is consistent 



with observations (e.g., Kormendy & Bender 2011). In addition, some stellar population in 
the outskirts (either on a disk or just at large radii of an elliptical galaxy) of AGN hosts 
may be unrelated to the preceding starburst and could be substantially different from bulge 



stars (e.g., Nolan et al. 2001). 
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Fig. 2. — shows the entire evolutionary process for an example merger of two gas-rich spirals 
of mass ~ 10^^ Mq each that eventually produces a powerlaw elliptical galaxy of velocity 
dispersion of 200km/s. This scenario is not limitd to merger events but encompasses any 
significant event triggering a starburst. Note that the time boundaries between difference 
consecutive phases are approximate and uncertain to within a factor of a few. The numbers 
in brown indicate the BH masses and the numbers in red indicate SFR. These numbers are 
very approximate and given mainly for illustration purpose. Clearly, given the complexity, 
one should expect large variations from case to case. 
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Fig. 3. — Top left panel: evolutionary growth tracks in the SMBH mass- Eddington ratio 
plane of an example SMBH of final mass 10^ with two cases of seed black mass of 10'' and 
10^ M0, respectively. A low mass cutoff for the IMF of 0.92 Mq that has a turnoff lifetime of 
10 Gyr is assumed. We assume that the SMBH accretion rate is proportional to the recycle 
gas return rate of the form oc t~^'^ Ciotti et al. (1991) capped at the Eddington rate with a 



radiative efficiency of e = 0.1, starting 200Myrs after the end of the starburst. Also indicated 
along each track are the times in Gyrs elapsed since the start of the accretion. Top right 
panel: the case for a low mass cutoff for the IMF of 1.4 Mq that has a turnoff lifetime of 
2 Gyr. Bottom panels: tracks for the cases in top panels but in the SMBH mass-luminosity 
plane. 
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7. Conclusions 

We have presented an alternative physical model that has the following characteristics 
for the coevolution of galaxy and SMBH. From the onset of a starburst to becoming a 
quiescent bulge (in the absence of any subsequent significant burst event) there are three 
distinct periods: 

(1) "Starburst Period": some significant event induces a starburst that probably lasts about 
10 — lO^yrs. The SMBH grows modestly during this period to possibly attain a mass that is 
up to order ten percent of its final mass. The feedback energy /momentum from the starburst, 
i.e., supernovae, drives the last patch of gas away and shuts down star formation. 

(2) "SMBH Prime Period": several hundred million years after the end of the starburst, aging 
low-to-intermediate mass stars, now in their post-main-sequence phases, start to return a 
substantial fraction of their stellar mass to the ISM. Because the rate of gas return from 
stars diminishes with time, the Eddington ratio of the SMBH decreases with time roughly 
as t~^-^. The SMBH growth is synchronous with star formation from recycled gas during 
this period. The accompanying star formation rate may also be substantial. The duration 
of this phase depends sensitively on the lower cutoff mass of the initial mass function (IMF) . 

(3) "Quiescent Bulge": on order of gigayear after the end of the starburst the elliptical 
galaxy is now truly red and dead - gas return rate is now negligible so both accretion to the 
central SMBH and residual star formation have ceased. It is possible that a disk may grow 
around the bulge later. The feeding of the central SMBH in the bulge of spiral galaxy during 
this period is not by overhead material from aging stars, rather by occasional objects that 
happen to be on some plunging orbits to be disrupted by the SMBH and form a short-lived 
accretion disc. Candidate objects may include molecular clouds or tidally disrupted stars. 

In this model, the end of starburst precedes the onset of prime SMBH growth by order 
of lOOMyr. Starburst is responsible for shutting down its own activities; AGN has little 
to do with it. AGN does provide self-regulation during its prime growth post-starburst 
period. An important feature of this model is that it is physically based and no significant 
fine tuning is required. The physical reason why the SMBH does not grow substantially 
in the starburst phase, although over-supplied with gas, is that only a very small central 
disc is gravitationally stable for gas accretion onto the SMBH, while all other regions are 
unstable and more conducive to star formation. The condition is just the opposite during 
the post-starburst phase where recycled gas dropout from aging stars returns slowly and 
can be more effectively accreted, so effective that self-regulation comes to play, energetically 
feasibly provided by the radiation pressure. 

Many comparisons between this physical model and extant observations are made and 
the model appears to be in very agreement with them, including the Mbh — Mbg relation. 
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This model predicts that the distribution of the Eddington ratio of AGNs in star-forming 
galaxies is lognormal, whereas that of AGNs in early type galaxies is a powerlaw, consistent 

. We predict that early-type galaxy hosts of 
high Eddingtion rate AGNs are expected to be light-blue to green in optical color, gradually 
evolving to the red sequences with decreasing AGN luminosity. 



with observations ( Kauffmann fc Heckman||2009 
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